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Abstract 

The one dimensional Hubbard model with nearest and (negative) 
next-nearest neighbour hopping has been studied with the density- 
matrix renormalization group (DMRG) method. A large region of 
ferromagnetism has been found for finite density and finite on-site 
interaction. 

PACS : 74.25.Ha, 75.10.Lp 

1 Introduction 

The Hubbard model was originally introduced to describe correlation effects 
in transition metals, such as, for example, the band ferromagnetism of Fe, 
Co and Ni. In mean field theory, one finds ferromagnetism whenever the 
Stoner criterion is satisfied, which, in the Hubbard model, leads to substantial 
regions of ferromagnetism in the phase diagram. When correlation effects are 
taken into account, however, a fully polarized ground state is quite difficult 
to find. For dimension d — 1, a theorem proven by Lieb and Mattis requires 
the ground state of the Hubbard model with near-neighbor hopping to be a 
singlet, completely excluding ferromagnetism H. For d > 2, ferromagnetism 
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is possible via the Nagaoka mechanism at very large interaction and near 
half-filling, but has not been found elsewhere in the phase diagram. Here 
we concentrate on the ID case which can be extended to circumvent the 
Lieb-Mattis theorem so that ferromagnetism is allowed. One way to do this 
is to add orbital degeneracy to the model in order to mimic the Hund's rule 
ferromagnetism found in atoms; another way is to add interaction terms such 
as nearest-neighbour Coulomb terms ||. Finally, one can change the band 
structure by introducing a "flat band" or, as is done here, by adding 
longer range hopping terms. 



2 Model and method 

We will study the ID Hubbard model with nearest and next-nearest neigh- 
bour hopping: 

L L 

H = - Yl (iicJ CT c i+ia + ^cL c i+2<7 + h.c.) + ^^n iT n a . (1) 
i=l,<r=U »=1 

Here c\ a creates an electron of spin a on site i, n ia = c\ a c i(7 , L is the system 
size, and U is the on-site Coulomb interaction. In the following we express 
all energies in units of t\ = 1 and only consider negative values of t%- Because 
a definite order of the particles is no longer enforced when t 2 ^ 0, the Lieb- 
Mattis theorem does not apply and, indeed, ferromagnetism has analytically 
been shown to exist at U = oo in three different limits: For one hole in a 
half-filled band Nagaoka ferromagnetism has been found Q; for j^l — * 0, it 
has been shown [H that the model is ferromagnetic for all densities; and for 
\t 2 \ > 0.25, where the band structure has two minima, Muller-Hartmann || 
has shown that the low density limit is ferromagnetic. These three limits are 
indicated in the schematic phase diagram shown in Fig. 1. In addition, for 
| £2 1 — ► 00 the model can be mapped onto two decoupled Hubbard chains, 
which cannot be ferromagnetic due to the Lieb-Mattis theorem. The aim of 
this work is to determine the extent of these domains of ferromagnetism in 
n and U, and whether or not these regimes are connected. 

Previous work using Lanczos exact diagonalization and variational tech- 
niques with various trial functions 0, |8| has shown that the ferromagnetic 
domain is large. However, the exact diagonalization calculations in this work 
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showed large finite size effects associated with closed shell effects in momen- 
tum space. Here we use the density matrix renormalization group (DMRG) 
method developed by White || and now widely used, to obtain results for 
much larger system sizes using open boundary conditions. 

3 Results 

Since the ground state energy obtained in the DMRG, Ed, is variational, it 
provides an upper bound for the exact ground state energy. We can analyse 
the stablity of the fully polarized state by comparing its energy, Ep, which 
is exactly known since the state has no double occupancy, with Ed- For a 
given value of t 2 and density n we can find a critical value U c below which the 
ferromagetic state is unstable. Since Ed can be found with high precision 
for this model, U c can be determined accurately. To check this we have 
calculated the total spin S in the DMRG ground state by evaluating 

(S 2 ) = E(S 4 S J ). (2) 

We find that the value of S goes smoothly from to the fully polarized 
value Smax & t the U c found by comparing the energies. It is difficult to 
determine whether the change in S is continous or not. The near-degeneracy 
of states with different S at the transition leads to mixing of states in the 
diagonalization step of the procedure, so that (S 2 ) no longer takes on definite 
discrete values. Using these two tools we have calculated the U — oo phase 
diagram shown in Fig. 2. The full dots are points at which the ground state 
is fully polarized (i.e. Ed > Ep and S ~ Smax ); the empty squares are 
points where the system is not magnetic (i.e. Ed < Ep and S ~ 0); and 
the dashes are points at which the ferromagnetism is not fully saturated (i.e. 
Ed > Ep and S ^ Smax). We can see that the previous analytical limits are 
well reproduced in this phase diagram and that the regions associated with 
the three mechanisms are connected. 

We next examine the behavior at finite U. In Fig. 3, we choose three 
representative values of t 2 (—0.1, —0.8, —2.0), and show U c as a function of the 
density n. At t 2 = —0.1, the band has only one minimum, at t 2 = —2.0 the 
ferromagnetic region does not occur for all densities at U = oo, and t 2 = —0.8 
is intermediate between the two regions. We see that for the two cases in 
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which there are two minima (1^1 > 0.25) in the band structure, there is a 
local minimum in U c at a given density n c . This density corresponds to the 
Fermi level being just at the top of the barrier between the two band minima, 
which leads to a high density of states, favourable for ferromagnetism. 

4 Conclusion 

We have found a ferromagnetic ground state in a large parameter regime in 
the ID Hubbard model with next-nearest neighbor hopping, showing that 
ferromagnetic regimes found in particular limits at U = oo extend to finite 
U and density. We also find that these different ferromagnetic regimes are 
connected. One interesting unresolved issue is whether the transition from 
a paramagnetic to a ferromagnetic state is continuous. This issue can be 
clarified in future work by looking at appropriate correlation functions in 
more detail. The present results confirm that the "critical density" n c in- 
troduced in previous work |7| does not represent a true phase boundary [§], 
but merely a density where the U value needed to stabilize ferromagnetism 
becomes small. 
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Figure captions 
Fig. 1 

Schematic U — oo ground state phase diagram in n — t 2 plane. 
Fig. 2 

The full dots show a ferromagetic groundstate, the open squares are for 
paramagnetic ground state, and the dashes are for points where doubts re- 
main. The energy and S where calculated on a system of size L = 30 at 
U = 10 6 . 

Fig. 3 

Critical value of U above which ferromagnetism is found for three different 
values of t 2 , namely t 2 = —0.2 (full dots), t 2 = —0.8 (open squares) and 
t 2 = —2 (stars). The lattice size is L = 50. 
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